Noise enhanced stability of periodically driven metastable states 
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We study the effect of noise enhanced stabihty of periodically driven metastable states in a system 
described by piecewise linear potential. We find that the growing of the average escape time with 
the intensity of the noise is depending on the initial condition of the system. We analytically obtain 
the condition for the noise enhanced stability effect and verify it by numerical simulations. 
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Escape from a metastable state is a phenomenon ob- 
served in several scientific areas. Among theni there are 
the theory of diffusion in sohds, chemical kinetics and 
transport in complex systems [y . 

The mean first passage time (MFPT) of a Brownian 
particle moving in potential fields usually decreases with 
noise intensity according to the Kramers formula H] or 
some universal scaling function of the system parameters 
[^,|| . The dependence on the noise intensity of the MFPT 
for metastable and unstable systems was revealed to have 
resonance character first noted by Hirsch et aLH and 
then observed in different physical systems P-|l 1[| . The 
most important conclusion of these studies is that the 
noise can modify the stability of the system in a counter- 
intuitive way. The system remains in the metastable 
state for a longer time than in the deterministic case and 
the escape time has a maximum at some noise intensity. 
Noise enhanced stability (NES) originally found numeri- 
cally by Dayan et al. H , was observed experimentally in 
a tunnel diode by Mantegna and Spagnolo |^. More re- 
cently it was found that the noise induced slowing down 
pO| and the noise-induced stabilization ||ll| are related 
to NES phenomenon g . 

Some questions arise from previous studies: (i) What 
is the reason of the increase of the average escape time 
with the noise intensity? (ii) What about the condition 
for which the NES effect takes place? To answer both 
questions we investigate the escape time from a period- 
ically driven metastable state for a piecewise linear po- 
tential. We study the nature of this phenomenon ana- 
lytically. We find that for fixed potential the decay time 
of unstable initial state can be dramatically increased by 
the presence of a small noise depending on the initial 
condition of the system. We obtain the condition for the 
NES effect, as an explicit relation between the driving 
frequency and the parameters of the potential. 

We consider the model of overdamped Brownian mo- 
tion described by the equation 
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dU{x) 
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where ^{t) is the white Gaussian noise with zero mean, 
{^{t)^{t + r)) = 2qS{T), F{t) is the dichotomous driving 
force and U{x) is a potential profile defined as 



U{x) 



CX), 


x^Q 


hx, 


Q <x<(. 


E -k{x- 


-i), £<x<b 



(2) 



with E = h£. Specifically we assume fc > and \F{t)\ < 
k. 

First we consider the system governed by Eq. (nh with 
F{t) — and potential (2) with arbitrary h. If ft. > 
{E > 0) the states at x < ^ become metastable. The 
exact expression of the MFPT from initial position xq to 
boundary b is known for the case F(t) — 



T{xo,b,q) = - / e 



u{z) 



-u(y) 
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where u{x) — U{x)/q is a dimentionless potential pro- 
file. Evidently in physical systems we can not observe 
the microscopic initial conditions. However the MFPT 
of Eq.(|) it is suflicicnt to obtain the MFPT with ar- 
bitrary macroscopic initial distribution by simple inte- 
gration. Therefore further we study T{xo,b,q) because 
it contains the full information about the system. If 
< xq < .^, the decay time r = ti for the potential 
profile of Eq. (0) is 
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If i? < xq < fe, the decay time r = r2 is 



T2{xa,b,q) = 
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(1) Here A = k{b - £), and AE = k{xo - £). 



These expressions show that at large noise intensity 
the decay time T{q) decreases with noise as 1/q for ar- 
bitrary h. When the noise intensity is small, q -^ \E\ 
the influence of the potential barrier becomes significant. 
For E < and h < 0, the barrier is absent, and the NES 
eff'ect, also known as noise delayed decay, appears when 
\h\ < fcandxo is near € Q. Indeed, for g < \E\,A, \E-A\ 
the decay time T2{i, b, q) grows with noise temperature q. 
When xq = 0, the decay time always decreases with q. 

In the case with potential barrier {E > 0, h > 0) the 
escape time depends on the initial position of the particle 
with respect to the potential barrier. When the particle 
is within the potential well {xq < £) the decay time of 
metastable state increases infinitely when g ^ 0, because 
if the noise is absent, the particle can never surmount the 
potential barrier. For q <^ E the decay time (^ coincides 
with the Kramers' time, which in this case reads 



Ti{xo,b,q) ~ Tk 



lih + fc) E/q 



(6) 



When i < Xq < b the initial state of the particle is un- 
stable. In the absence of noise, the escape time from this 
unstable state is a finite value: T2(xo, b, 0) — (b — Xo)/k, 
which does not depend on the potential well. When we 
add the noise the influence of potential well becomes im- 
portant. It follows from the exact expression of Eq.(P) 
that the MFPT rises to infinity ii AE < E 
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while for the case £,{t) — we have the decay time ob- 
tained from the deterministic Eq. (P, i. e. the MFPT has 
a singularity at q = 0, when AE < E. From a physical 
point of view this singularity can be explained as fol- 
lows: When the particle is initially located in the region 
£ < xo < b, a. small quantity of noise added in the sys- 
tem can eventually push the particle into potential well. 
Then, the particle will be trapped there for a long time 
because the well is very deep. This type of trajectories of 
the Brownian particles therefore leads to a big "tail" in 
the First Passage Time Distribution (FPTD) w(t). If the 
potential well is very deep, namely, E > AE the trapping 
time is so long that the integral for MFPT 



r = / tw{t)dt 
/o 
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diverges when q ^ {). The FPTD obeys the backward 
Fokker-Planck equation and can be obtained for piece- 
wise linear potential using the Laplace transform method 
[ p2[ . The Laplace transform of the FPTD for our poten- 
tial, when h = k and £ < xq < b, reads 
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where 
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fi = p + c, (3 — k£/q. Using the limit theorems of Laplace 
transform, we can obtain from Eq. (S) the asymptotic ex- 
pression for i — > oo and q <^ E 
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where e{x) = (x - 2^)/fc, to = 9{b) ~ e{xo), h = 
9{b) + 9{xo), and Tk — Aq/k^. The time dependence of 
Was{t) is shown in Fig. 1 for three different values of the 
initial position xq = 2.5 {AE > E), xq = 2 {AE = E), 
and xq — 1.2 {AE < E). One can see that the tail of the 
FPTD rises when AE decreases. If AE > E, the trap- 
ping time in the well is not very long and the integral of 
Eq.(||) always converges. Nevertheless, the average decay 
time increases with small noise, reach the maximum and, 
then, decreases. The plots of T2{q) for various relations 
between AE and E are shown in Fig. 2. The main con- 
clusion from the above analysis is that the strong effect 
of NES can appear for the fixed potential profile with 
barrier, if the initial probability distribution is located 
within the interval {I, 6), i. e. in an unstable state beyond 
the potential well. The physical system can be brought in 
this nonequilibrium state by a sudden change of control 
parameters. Examples of such situations include spin- 
odal decomposition in the dynamics of phase transitions 
and the process of laser switch-on [Hjl^ . Such relaxation 
processes in the systems which are far from equilibrium 
attract now a great deal of attention ||l^,^,Q . 

The main aim of this paper however is to study the 
escape from the metastable state with an initial distri- 
bution located within the potential well in the presence 
of periodical driving. Further we will apply the above 
results, obtained in the static case, for analysis of pe- 
riodical force effect. Let's consider the same potential 
profile U{x) of Eq.(g) but with h = 0. The driving force 
is F{t) = ai'{t), where v{t) is the dichotomous signal 
switching between ±1 with period T and a is the am- 
plitude. We choose v{t) = —1 (i. e. the barrier is ab- 
sent) for the first half of period. The exact expression for 
the MFPT when potential varies with time is unknown. 
Three recent papers |lj] develop a theory of escape rates 
for periodically driven systems. Smelyanskiy et al. con- 
sider the case of high potential barrier and small driving 
amplitude. Lehmann et al. consider the regime of strong 
and moderately fast driving. Maier and Stein analyse 
the crossover regime. In all these cases the deterministic 
escape time is infinity. In the present work we consider 
the different regime of strong and moderately slow mod- 
ulation when the deterministic escape time is finite. 



Therefore we first start our investigation from the de- 
terministic case (7 = and secondly we define the con- 
dition for NES effect when noise intensity is small. If 
xq — 0, the Eq. (||) has a periodic solution in the deter- 
ministic regime for T < 2£/a. In this case x{t) < i for 
any t and the particle always remains in the metastable 
state. If the period is 



e ab + ki 
2- <T <2- -. 

a a(a + k) 



T > 2£/a 



(10) 



the particle surmounts the potential barrier at time 
t = i/a. To obtain NES effect we should consider only 
the cases when the states are unstable without noise [||. 
Consequently the first condition for the NES is given by 
Eq.(no[). It follows from the above analysis that the de- 
cay time in the presence of noise strongly depends on po- 
tential barrier, namely, the barrier is responsible for the 
strong increasing of decay time with noise. The exact ex- 
pressions for fixed potential show that increasing of decay 
time in the case without barrier is much smaller and it 
appears only if the particle is near the point xq — i (see 
Eq.s (|j) and ^)). Therefore it is important to consider 
two cases: (i) t(0, 6, 0) < T/2 and (ii) t(0, 6, 0) > T/2. 
In the first case the modulation frequency is so low that 
during the entire process of decay the potential profile 
^{x,t) has no barrier. In this case the average escape 
time is not increased by the noise (see Eq. (Q) at xq = 
and h — —a). In case of t(0, b, 0) > T/2, average escape 
time can be represented as follows 



r(0,6,0) = |+T(a;(r/2),6,0) 



(11) 



where x{T/2) is the position of the particle at time t ~ 
T/2. The conditions <^) and r(0,6,0) > T/2 together 
mean that 



I < x(T/2) < b. 



(12) 



Now we add a small quantity of noise into the system. 
The MFPT r(0, b, q) can be written as 



r(0, b, q) = t(0, x^,q) + T{x^,b, q), 



(13) 



where Xi is an arbitrary point between and b. For 
Xi = x{T/2) and very small noise, the first term in the 
right hand side of Eq. ([l^ ) is approximately equal to the 
deterministic time: t(0,x(T/2), g) ~ t(0, a;(r/2), 0) = 
T/2, because the MFPT varies smoothly with noise, 
when the barrier is absent (See Eq. (Q)). In this case: 
r(0,6,(7) ~ T/2 + t(x(T/2),6,9), and r(x(T/2), 6, 9) » 
t(x(T/2), 6, 0) because of the potential barrier, which 
makes the average escape time very large just for q —> 
(see Eq.(R)). As a result the decay time t(0, 6, (7) will 
increase with q and the NES appears. (The decay time 
r(0, 6, q) will not grow infinitely at q — > O because the 
barrier exists only a half of period.) Thus, we may con- 
clude that the inequality (O) must be the condition for 
NES. This inequality can be rewritten as follows: 



(14) 



The inequality (^^ and the condition a < k give the 
area on the (T, a) plane where the NES effect takes place. 
In Fig. 3 we show this area for fc = 1, ^ = 2, and b — 7, 
and the results of numerical simulations (shaded area). 
We perform 3000 different realizations of the decay pro- 
cess x{t) to determine the average escape time for each 
couple of values of the amplitude a and the period T of 
the driving force. We consider more than 100 points on 
the (T, a) plane. We find that within the area defined 
by inequality ( |l4| ) the NES effect is very strong: the av- 
erage escape time increases more than 10% above the 
deterministic escape time. Outside this area and below 
the lower boundary the deterministic decay time becomes 
infinite and the NES disappears. In the presence of noise 
we obtain Kramers-like behaviour. This case was stud- 
ied in detail by Lehmann et al. |l4|. Inside the area the 
magnitude of the NES effect decreases from the lower 
to the upper boundary. Above the upper boundary the 
NES effect decreases sharply. When the period T and the 
amplitude a of the driving force are chosen near the up- 
per boundary of (M), the potential barrier is very small 
or absent during the process of decay. It explain why 
the effect is very small, when we are near this boundary. 
We also carried out simulations for a > k and found the 
NES effect. This parameter region however gives less in- 
formation about the mechanisms of the NES effect from 
viewpoint of the interplay between the regular, random, 
and periodical forces. In this region in fact the determin- 
istic motion of the particle is characterized by oscillations 
and the driving force prevails the regular one described 
by the potential U{x). 

Therefore we conclude that the numerical simulations 
are in good agreement with the theory. The main mech- 
anism of the NES is defined correctly: it is the role of 
the potential barrier which appears after the particle has 
crossed the point £ of maximal potential. Consequently 
the inequality dl2) is the most general condition for NES 
effect because it can be applied to a system described by 
an arbitrary potential with metastable state and where 
£ is the x-coordinate of the maximum of the barrier ap- 
pearing a,t t — T/2. The mechanism of NES explains 
why the EFT distributions obtained in simulations and 
in experiments are multipeaked, periodic and with an 
exponential time decaying envelop P,p|. The peaks ap- 
pear only for small noise intensity, where the NES effect 
occurs. The first peak corresponds to the deterministic 
escape time. The second peak arises because the small 
noise provides the above considered inverse probability 
current which moves some particles into the potential 
well. The returned particles can escape only in one pe- 
riod. Therefore the second peak is one period apart from 
the first one. After each period we have the same physi- 
cal situation and as a consequence fewer particles go back 



into the potential well. Therefore the probability peaks 
have period T and they decrease with time. The proba- 
bilities of escape are independent and equal for successive 
oscillations of the potential. So if the escape probability 
per oscillation is p, the probability to escape at the n-th 
cycle is (1 —p)"^^p ~ pe^e^^r where —a = ln(l —p) and 
n ~ t/T. Therefore the magnitude of the FPT peaks are 
exponentially decreasing with time. 
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E dot-dashed line), xo ^ 2 {AE = E full line), and xo = 1.2 
{AE < E dashed line). The parameters are: b — 3,k — h = 
1,1 = l,g = 0.01. 

FIG. 2. The normalized decay time T2{q)/Td (with Td the 
deterministic time) versus q/E for the same parameters of 
Fig.l. Inset: the potential U{x) of Eq. (2). 
FIG. 3. The shaded area is the region of the plane (InT, a) 
where the NES effect is very strong: the average escape time 
is greater than 10% above the deterministic escape time. The 
parameters are: b = 7,k = 1,1 = 2. Inset: the average escape 
time versus the noise intensity for a = 0.3 and T — 13.5. The 
dashed line indicates the deterministic escape time. 
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FIG. 1. Asymptotic behaviour of the FPTD Was{t) versus 
time t for three values of the initial position: Xo = 2.5 (A_E > 
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